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In this paper I study the posibility of inducing a single-electron current by rotating a non-magnetic conducting 
rod with a small tunnel junction immerse in a uniform magnetic field perpendicular to the plane of motion. I 
show first, by using a thermodynamic approach, the conditions needed to pump electrons around the mechanical 
device in the Coulomb blockade regime. I then use a density matrix approach to describe the dynamics of the 
single-charge transport including many-body effects. The theory shows that it is possible to have single-electron 
tunneling (SET) oscillations at low temperatures by satisfying conditions similar to the Coulomb blockade 
systems. 
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INTRODUCTION 

The field of single-electronics started when new effects due to the quantization of charge in ultrasmall tunnel junctions, both 
in the superconducting and the normal state, where predicted by Averin and Likharev [1 1. The theory of Averin and Likharev 
considers a tunnel junction which is biased by an externally fixed current i and whose voltage V is measured by a very high 
impedance voltmeter with metallic shunt conductance Gs- A tunnel junction consists of two conducting electrodes separated 
by a thin layer of insulating material and is characterized by its capacitance C and tunnel resistance Rt- When a voltage is 
applied to the small capacitance tunnel junction, the charge will flow continuously through the conductor and it will accumulate 
on the surface of the electrode against the insulating layer of the junction (the adjacent electrode will have equal but opposite 
surface charge). On the other hand, the insulating layer is thin enough for electrons to tunnel through. The state of the junction 
is described by the surface charge Q (which is a continuous variable) and the electrons n that have tunnel through the insulating 
layer (which is a discrete variable). Averin and Likharev predicted that if the charge Q at the junction is greater than |e|/2, an 
electron can tunnel through the junction in a particular direction, subtracting |e| from Q. Likewise, if Q is less than — |e|/2, 
an electron can tunnel through the junction in opposite direction, adding |e| to Q. But if Q is less than |e|/2 and greater than 
— |e|/2, tunneling in any direction would increase the energy of the system, hence tunneling will not occur. This suppression 
of tunneling is known today as the Coulomb blockade Q. The physical origin of the Coulomb blockade of single-electron 
tunneling (SET) is quite simple. In a current-biased junction, each tunneling event leads to a change of the electrostatic energy 
of the system given by AE = e(Q±e/2)/C. If the initial charge Q is within the range — | e | /2 < Q < |e|/2, the energy change 
AE is positive and at low temperatures tunneling events are impossible. On the other hand, if |Q| > |e|/2, tunneling is possible 
because this process reduces the electrostatic energy. An interesting prediction of Averin and Likharev was the SET oscillations 
in the voltage across the junction Q]|2]. Due to the Coulomb blockade of tunneling the charge Q on the junction accumulates 
until its threshold value e/2 is reached and then the junction is recharged by the externally fixed current. This whole process 
repeats itself with a frequency v = i/e |3l. 

The purpose of this article is to show that there are SET oscillations without an external applied current in a small tunnel junction. 
In this case, the SET oscillations are driven by the Lorentz force due to the rotation of a conductor with a small tunnel junction 
and an applied external magnetic field. In addition, this mechanically driven device is proposed as a transducer of motion into 
electricity. The paper is organized as follows. First I will start by giving a thermodynamic formulation of the problem and the 
basic relations of the theory. Then I will analyze the system using a density matrix approach to include many-body effects and 
show the SET oscillations in the system. The conclusions are summarized in the last section. 

THERMODYNAMIC FORMULATION 

Consider the Faraday generator shown in Fig. ([TJ, where a conducting rod of length £ rotates with constant angular velocity 
uj in a constant magnetic field that is perpendicular to the plane of motion. The rod completes the circuit, with one contact 
point on one end of the rod and the other on the circular rim. The circuit containing the galvanometer is completed by an open 
wire structure with a switch. For the case when the switch is open there will be no electromotive force, however we know from 
elementary electrodynamic courses that charge will pile up at the two ends of the rod and will produce an electric field that 
balances the Lorentz force felt by the moving charges inside the conductor (See Fig. (|T]i(a)). 

The free energy of the system is given by T = E — fl ■ B, where E is the total energy of the rotating rod, fl is the magnetic 
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dipole moment of the circulating charge at the end of the rod and B is the constant magnetic field perpendicular to the plane of 
motion. Taking the magnetic field as B = —Bz, the free energy is given by 



T — Ei. 
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where Ei„ is the internal energy of the rod, / denotes the moment of inertia for the rod with respect to the axis of rotation and 
we have taken the current as i = Qu/2ir. It should be remembered that rotation in general changes the distribution of mass in 
the body, and so the moment of inertia and internal energy of the body are in general functions of u> (4) . 

When the switch is turned on there is an electromotive force £ in the circuit and hence charge will be circulating around the 
circuit. The free energy of the system in this case is given by 



T — Ei, 
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where AQ is the amount of charge circulating around the circuit and £AQ is the work done by the system. The change in free 
energy is obtained by subtracting Eq.Q from Eq.Q, which gives us 
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if the system is in thermodynamic equilibrium then AT 7 = and we obtain £ = —Buj£ 2 /2, as we know from elementary 
electrodynamic courses 0. Note that £ < 0. 

Now I will consider the case when there is a tunnel junction at position r with thickness Sr in the conducting rod as shown in 
Fig.Q. For this case, even if the switch is turned on there will be no current flowing through the circuit because the charge 
Q will accumulate on the surface of the electrode against the insulating layer as depicted in Fig. ([2]). Nevertheless, quantum 
mechanically speaking there is a probability for the charge to tunnel through the junction. The free energy of the system before 
quantum tunneling is given by T = E — fl ■ B — fa ■ B — fa ■ B, where /2 and /Ji(2) corresponds to the magnetic dipole moment 
of the charge accumulated at the end of the rod and against both sides of the insulating layer, and E = Iuj 2 /2 + Q 2 /2C, where 
C is the capacitance of the tunnel junction. Therefore, the free energy of the system before quantum tunneling is given by 
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Eq. Q can be written in the following form 
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FIG. 2: Schematic diagram for the rotation of the conducting rod with a tunnel junction of thickness Sr and a uniform magnetic field pointing 
into the page. Note how the charge accumulates on the surface of the electrode against the insulating layer. For this case current will only flow 
when a tunnel event occurs. 



When there is quantum tunneling there is a change in charge by ±|e| and an electromotive force £ in the circuit which causes a 
change in the free energy given by 

V 2 (Q'±|e|) 2 

F= 2 U)+ 2C +g|e| " (?) 
The change in free energy is obtained by subtracting Eq.|7]) from Eq.Q, which gives us 



2C V |e| , 

Assuming £ < from our previous result, we see from equation ([8]l that a tunnel event becomes energetically favorable and a 
current i flows throughout the circuit when Q' > When a tunnel event takes place the charge Q' will change by — |e| and 
after a time \e\/i the rotational motion of the conducting rod immerse in the constant magnetic field will recharged the junction 
and another tunnel event will take place. As a result the tunneling events will occur periodically with frequency v = i/\e\. 
For the particular case in which the tunnel junction of thickness S lies exactly in the middle of the conducting rod, i.e. r = 
1/2 — 8/2 and r + Sr = i/2 + 5/2, then a tunnel event becomes energetically favorable when 

Eq. |9]l can be expressed in terms of the electrostatic voltage in the following way 

lei Bu£(£-6) 

V> 2CJ + ^2-^- (10) 



Since the maximum voltage allowed for the system is given by Bu>£ 2 /2, then Eq. ( 10 1 reads 

2 2C 2 K ' 

Eq. <[□) gives us the following restriction Buj£S > \e\/C. Using typical values of the tunnel junction capacitance C « 
3 x 10~ i5 F and tunnel thickness S = 10A J2j , we need BloI > 10 5 V/m, to satisfy the restriction condition. If we have a 
magnetic field of B = IT and £ — 1cm, then we will need a rotational frequency of around v r > 10MHz. The current 
delivered by the single-electron Faraday generator for this rotational frequency is around i « IpA 



MANY-BODY EFFECTS 



To study in more detail the dynamics of the charge transport we need the total Hamiltonian of the system depicted in Fig. ©, 
which is given by: 

U = F{Q')+H T + H 1 +H 2 +H s -i<$>. (12) 
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The first term in Eq. ( 12 1 represents the free energy of the system which is given by Eq. (^J. The charge operator Q' can be 
expressed via Fermion operators 



\ fcl fc 2 / 



where cl and c k are the electron creation and annihilation operators and Qo = Gui 
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is a constant term. The 



second term in Eq. ( 12 1 represents the tunneling Hamiltonian which is given by 
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where the summation is carried out over all states k within the electrodes 1 and 2 and T kl k 2 is the tunneling rate across the 
junction. The Hamiltonians Hi, H% and Hs describe the energy of the internal degrees of freedom ki, and k$ of the two 



electrodes of the junction and of the shunt Gs, respectively. The last term in Eq. ( 12 1 is the operator of the magnetic flux defined 
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where £ is the electromotive force (emf) around the circuit. Note that Eq.( 12 1 corresponds exactly to the basic Hamiltonian 
given by Averin and Likharev. The only new feature is the shift of the charge operator Q' = Q — Qo arising from the magnetic 
dipole interaction between the spinning charge and the external magnetic field. 

Restricting ourselves to the case when the current through the junction and shunt are not too large we can consider them as 
perturbations and one can write an explicit time evolution equation for the density matrix to describe the junction properties. 
Following the pioneering work of Averin and Likharev [ 1 ], and assuming that Gs, Gt « 4e 2 //i, the resulting master equation 
is given by 
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where f(Q', t) is the classical probability distribution and Ft and F$ are contributions due to the tunneling and shunt current, 
respectively, and are given by 



F T (Q') = r+(Q'-e)f(Q'-e,t)+r-(Q' + e)f(Q , + e,t) 
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where are the tunneling rates for forward (plus sign) and backward (minus sign) single electron tunneling over the junction 
and can be expressed as 



r ± (Q / ) = -t(AF ± /e) 
e 



where i is the d.c. induced current and 



AJ"± =±g(e/2±Q'). 



(19) 
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Note that the master equation given in Eq. ( [T6| > corresponds exactly to the master equation given in Averin and Likharev paper 
1 1 1, the only difference is that there is no external current or bias. 



If one looks at the regime — e/2 < Q' < e/2, the tunneling contributions may be neglected, i.e. Ft = 0, and Eq.( 16 1 can be 
expressed as 
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where V(t) is the time dependent voltage across the tunnel junction and is given by |6| 
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Making the substitution F(Q', t) = f(Q', t) + V(t)/ k B T, we end up with a reaction-diffusion equation given by 

OF d 2 F 1 dV 

^t- GskBT dQ^ + k^rHt- 



For the system with constant electrostatic potential Vq, i.e. no tunneling, the solution to Eq. (23 1 is (TJ 
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where the first term in Eq. (24 1 is the solution to the master equation for f(Q', t), which properly normalized can be expressed 
as 
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Eq. (25 i represents a Gaussian probability packet describing the distribution of charge Q that will move due to the Lorentz 
force until Q' > e/2, at this point the rate r~(Q') becomes nonvanishing and this leads to a rapid decay of the packet, i.e. 
f(Q', t — t-r) — 0, where — C/Gt is the time when a tunneling event occurs. In this regime the tunneling event leads to a 
noticeable change in the voltage across the junction, i.e. AV = ±e/C, which is described in Eq. (23 1 by the last term, i.e. 
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The general solution to Eq. ( 26 1 is just a shift of the solution given in Eq. ( 25 1 to the starting time t = tx, i.e. 

f(n , _ e(t-t T ) -Q' 2 
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where 8(i — tx) is the Heaviside funtion. It is evident from Eq.(27i that the whole process of the Gaussian packet formation 
repeats periodically with every tunneling event showing the periodic SET oscillations in the voltage across the junction. 



CONCLUSIONS 



The main contribution of this article is to show that there can be SET oscillations across a tunnel junction without an externally 
applied current source. This result is in contrast to the system analyzed by Averin and Likharev where an external fixed current is 
always present. The thermodynamic and microscopic derivation shows how a single-electron current can be induced by rotating 
a conducting rod with a small tunnel junction in the presence of a uniform magnetic field perpendicular to the plane of motion. 
An estimate of the current delivered by the single-electron Faraday generator for rotational frequencies of v r ps 10MHz is 
i « lpA. Thus, this device could serve as a fundamental standard of d.c. current. 
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